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METHOD AND SYSTEM FOR PROCESSING VASCULAR RADIOGRAPHIC 
IMAGES WHICH HAVE BEEN RECONSTRUCTED BY THREE-DIMENSIONAL 
MODELLING 

CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] The entire contents of French Patent Application No. 0100737 
filed January 1 9, 2001 are hereby incorporated by reference. 

BACKGROUND OF THE INVENTION 

** 
O 

[0002] The present invention is directed to a method and system for 
I** processing radiographic images which have been reconstructed by three-dimensional 
]T modelling. In particular, the present invention is applicable to vascular images which 
'1 '' includes an implanted element, such as a prosthesis. 

h 

w 

[0003] There are known numerous three-dimensional imaging methods 
p and systems which take two-dimensional images obtained, for example, by X-ray 
r '* fluoroscopy, to produce 3D models of an object, i.e., a patient that it is desired to be 
observe. In particular, there are known 3D angiography systems which, by X-ray 
fluoroscopy, reconstruct 3D models of vessels on which a procedure is to be carried out, 
for example in order to treat arterial stenosis. In 3D angiography, there are three 
complementary types of image likely to be obtained, namely: 

[0004] a reconstructed model known as a "subtracted" model, indicating 
the vascular elements ("lumens") alone, without the calcified elements and the 
endovascular prostheses; 



[0005] a reconstructed model known as a "masked" model, identifying 
the calcified elements and the prostheses, but not the vascular elements; and 
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[0006] a reconstructed mode] known as an "opacified" model, 
identifying the vascular elements, the calcified elements and the prostheses, but without 
it being very easy to distinguish the various elements, the image obtained being 
relatively difficult to interpret. 

[0007] Techniques for producing these three types of modelled images, 
are described in French Patent Application No. 001 1486. 

BRIEF DESCRIPTION OF THE INVENTION 

[0008] An embodiment of the invention is system and method for 
processing a vascular radiography image which has been reconstructed by three- 
dimensional modelling comprising: 

[0009] (a) determine a three-dimensional model to be known as the 
masked model which features the calcified elements and an implanted element, but not 
vascular element; 

[0010] (b) determine a three-dimensional model to be known as the 
subtracted model, which features the vascular element alone; 

[001 1] (c) merging the two models and weighting their voxels so as to 
increase the contrast between the image of the masked model and the image of the 
subtracted model; and 

[0012] (d) summing the voxels thus weighted. 
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BRIEF DESCRIPTION OF THE DRAWING 

[0013] An embodiment of the invention will become further apparent 
from the description which follows, which is purely illustrative and non-limiting and is 
to be read with reference to the appended drawings in which: 

[0014] Figure 1 schematically depicts a system for implementing a 
method according to an embodiment of the invention; 

[0015] Figure 2 is a flow diagram schematically illustrating the various 
steps in the processing in an embodiment of the invention; 

[0016] Figure 3 is an example of a possible curve for one of the 
functions used in the processing; 

[0017] Figures 4a and 4b illustrate a 3D model and a sectional view 
obtained from it, in the case of implementation of a method of one embodiment of the 
invention; and 

[001 8] Figures 5a and 5b illustrate another exemplary embodiment. 

DETAILED DESCRIPTION OF THE INVENTION 

[0019] In this embodiment, it is assumed that there is available a set of 
two-dimensional angiography images obtained around a given anatomical region of a 
patient and from these images it is possible to reconstruct three-dimensional models of 
the anatomical region and, in particular, subtracted and masked three-dimensional 
models. These angiography images may, for example, be obtained by x-ray fluoroscopy, 
etc. The images are stored and processed in a processing unit 5, which is connected to 
interface means 6 which, in particular, allow the radiography images to be displayed 
(Figure 1). 
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[0020] In an embodiment of the processing may be carried out in four 
successive steps. 

[0021] In a first step (step 1 in Figure 2), the user defines a volume 
which constitutes the region of interest. The region of interest will, in particular, 
preferably be defined so that it contains only a limited number of bones while at the 
same time fully including the portion or portions of blood vessel the user wishes to 
view, with any calcified elements or prosthesis they may have. In particular, the user 
may select, from a displayed 3D model the vessels to view and for the system to 
automatically determine the limits of these vessels. As an alternative, the user can make 
a selection by magnifying the delimited regions. 

S 

[0022] In a second step (step 2), the system automatically reconstructs 
f,k the three-dimensional models which correspond to the subtracted reconstruction and the 
I- masked reconstruction for the selected region. 

BT 

[0023] In a third step (step 3), the system performs a 3D merger 
P processing operation on the two models thus obtained. 

PJ [0024] And finally, in a fourth step (step 4), the merged 3D image thus 

obtained is displayed on the interface means 6. 

[0025] The merger of the two models is substantially meaningless 
unless the subtracted and masked 3D models have been obtained from the same 
angiographic acquisition and correspond to the same 3D frame of reference. 

[0026] The merger processing operation is performed voxel by voxel, 
calculating the intensity (attenuation coefficient, Mer(v) of a voxel v of the merged 
image using the formula: 
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[0027] Mer(v) = V 0 . Sub(v) +a.Y b [Mask(v)] 
Mean(Sub) 

[0028] where: 

[0029] v represents the coordinates of the voxel; 

[0030] Sub(v) is the intensity (that is to say the attenuation coefficient) 
of the voxel v of the subtracted model; 

[003 1 ] Mask(v) is the intensity (that is to say the attenuation coefficient) 
of the voxel v of the masked 3D model; 

[0032] Mean(Sub) is the mean intensity calculated either over the entire 
volume considered, or vessel by vessel - using automatic determination of the limits of 
the vessels or vessel portion per vessel portion, or alternatively still, along straight-line 
portions which constitute the main directions of the vessel; 

[0033] V 0 is a predetermined constant which represents the "desired" 
mean intensity for the depiction of the vessel in the merged image; 

[0034] a is a constant of predetermined value which has a value higher 
than 1 (for example 3 0) so as to increase the contrast between the structures present in 
the unsubtracted model (vascular prosthesis, calcification, etc.) with respect to the blood 
vessels; 

[0035] Y b is a one-dimensional function which is, for example, a linear 
straight line, but which may be a more complicated function (c.f. Figure 3) which 
depends on a threshold value b and comprises three parts: 



[0036] 



Yb is zero in the interval [0,b]; 
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[0037] Y b varies linearly between b and a value b+b 0 , the 

linearity coefficient being equal to 1 ; 

[0038] Y b is, in the interval [b+b 0 ; +«], a function which 

increases less markedly than the linear function of intensity used between b and b+b 0 ; 
this function is for example the V function adjusted to ensure continuity of the first-order 
derivative of the function Y b , this third part making it possible to avoid the effects of 
saturation when highly absorbent metal elements are present in the masked image. 

[0039] The role of the threshold b is to eliminate the ambient noise 
present in the reconstructed or subtracted model before increasing its contrast. 

-g • 

% [0040] One possible robust automatic estimate of this threshold is 

m 

y= obtained as follows: 

#J 

[0041] calculating a smoothed histogram of the masked model, 

O and doing so over the selected region for treatment, smoothing being achieved by taking 
p a mean over a smoothing window of a predetermined size and, 

53 

M [0042] determining the minimum value of the first-order 

derivative of this histogram, the intensity thus obtained being the threshold b. 

[0043] Multiplying by the coefficient "a" is advantageously not 
performed until after filtering has been performed using the threshold b. 

[0044] Examples of merged images obtained in the way which has just 
been described are given in Figures 4a and 4b, and in 5a and 5b. The images in Figures 
4a and 4b being, respectively, views in section and perspective revealing calcifications 
around the vesselsl. The images of Figures 5a and 5b being, respectively, views in 
section and perspective revealing vascular prosthesis (stents). Merged three- 
dimensional images in which the calcifications and the protheses show up particularly 
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clearly are thus available. The calcifications or the stents show up particularly clearly in 
these images, making them easier to interpret. 

[0045] Of course, it is possible for the user, during the procedure, to 
alter certain parameters and in particular to modify the parameters b and a to adjust the 
contrasts of the various merged parts and control the brightness of the unsubtracted 
images. 



[0046] Images thus obtained are similar to those that can be obtained by 
CT angiography. A two-colour depiction may also be envisaged. 

%l [0047] Various modifications in structure and/or steps and/or function 

O may be made by one skilled in the art without departing form the scope and extent of the 

■ 

y s invention as recited in the claims. 
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